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A single crystal of the Co^^—based pyrochlore NaCaCo 2 F 7 was studied by inelastic neutron scat¬ 
tering. This frustrated magnet with quenched exchange disorder remains in a strongly correlated 
paramagnetic state down to one 60th of the Curie-Weiss temperature. Below T/ = 2.4 K, diffuse 
elastic scattering develops and comprises 30 ± 10% of the total magnetic scattering, as expected for 
e/eff = 1/2 moments frozen on a time scale that exceeds h/SE=3.S ps. The diffuse scattering is con¬ 
sistent with short range XY antiferromagnetism with a correlation length of 16 A. The momentum 
(Q) dependence of the inelastic intensity indicates relaxing AT-like antiferromagnetic clusters at en¬ 
ergies below ~ 5.5 meV, and collinear antiferromagnetic fluctuations above this energy. The relevant 
XY conhgurations form a continuous manifold of symmetry-related states. Contrary to well-known 
models that produce this continuous manifold, order-by-disorder does not select an ordered state 
in NaCaCo 2 F 7 despite evidence for weak (~ 12%) exchange disorder. Instead, NaCaCo 2 F 7 freezes 
into short range ordered clusters that span this manifold. 

I. INTRODUCTION 


FIG. 1. a) Elastic magnetic scattering from NaCaCo 2 F 7 ob¬ 
tained by subtracting elastic scattering at T =14 K from 
that at 1.7 K. Two diffuse components are visible: broadened 
Bragg peaks occurring at nuclear-allowed positions, and an 
extended “zig-zag” component that is also apparent at finite 
energy transfers [Fig. 4 b)]. The data were acquired at fixed 
monitor counts; we report approximate counting time per set¬ 
ting (10 seconds). The total time required to collect the data 
shown here was 32 hours, with 16 hours per temperature, b) 
The two basis states of the F 5 IR {XY states), combinations 
of which form the short range order in NaCaCo 2 F 7 . 


while quenched disorder is predicted to favor a coplanar 
LRO state in the same model.The role of quenched 
disorder for the HAFM pyrochlore model has been stud¬ 
ied in the past by including a distribution of exchange 
interactions spanning J ± A in the HAFM Hamiltonian, 
H = JijSi • Sj. In the limit of weak disorder, 
A << J, where J is the mean exchange interaction, the 
spins are expected to form locally collinear antiferromag¬ 
netic correlations and the system eventually freezes at a 
temperature Tf ^ 


The spin liquid state of the Heisenberg antiferromag- 
net (HAFM) on the pyrochlore lattice supports fluc¬ 
tuations within an extensively degenerate ground state 
manifold consisting of correlated, yet disordered, spin 
configurations.^'^ This beautiful state of matter arises 
from a perfect frustration of antiferromagnetic (AFM) 
interactions on the corner sharing tetrahedra that com¬ 
prise the pyrochlore lattice. However, the spin liquid 
is extremely susceptible to small perturbations that can 
reduce the ground state degeneracy and lower the free 
energy. The manner in which the spin liquid is modified 
in real materials with deviations from ideal Heisenberg 
exchange is thus a rich field of study, with many possi¬ 
ble outcomes depending on the relevant perturbations.^ 
In particular, the role of fluctuations in selecting sub¬ 
sets of the ground state manifold must often be consid¬ 
ered. Thermal and quantum fluctuations that are softer 
for certain spin configurations can, in some cases, select 
long range ordered (LRO) states in a mechanism called 
order-by-disorder. 

Quenched disorder, in the form of vacancies or bond 
disorder (i.e., local variations in the strength of the spin- 
spin interactions), also produces order-by-disorder, as de¬ 
scribed in the pioneering work by Villain^^ and later 
studied in detail by others.Quenched disorder 
can compete with thermal fluctuations to determine the 
ordered state. An important recent example is the XY 
antiferromagnetic pyrochlore material Er2Ti2 07. For the 
pseudospin ^ model believed to be appropriate for this 
material, thermal and quantum order-by-disorder have 
been shown to select a non-coplanar LRO state, 
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Here we report on the nature of static and dynamic 
spin correlations in the recently discovered pyrochlore 
material, NaCaCo 2 F 7 , which has been synthesized in sin¬ 
gle crystal form via the optical floating zone method. 
The inherent local disorder arising from the mixed-charge 
H-site (Na+/Ca^+) is expected to yield exchange dis¬ 
order. The availability of large single crystals has al¬ 
lowed us to measure the full dynamic structure factor for 
NaCaCo 2 F 7 in the high symmetry [HHL] plane. Our 
measurements indicate NaCaCo 2 F 7 adheres to the expec¬ 
tations for the HAFM with weak exchange disorder in the 
high energy limit {E > 5.5 meV), but at low energies it 
displays frozen short range correlations and relaxational 
dynamics associated with an easy plane (XT-like) man¬ 
ifold. 

In NaCaCo 2 F 7 (space group Fd3m, a =10.4056(2) A 
at T = 295 K^^) the H-site of the pyrochlore lattice is 
occupied by Na+ and Ca^+ ions with equal concentration 
in a disordered configuration. Thermodynamic magnetic 
properties of NaCaCo 2 F 7 evidence a spin-freezing tran¬ 
sition at Tf ^ 2.4 The low freezing temperature 
indicates weak exchange disorder with A/J ^ 0.12 (as¬ 
suming that Tf ^ A, as in the models of Refs. 16 and 
17 and J ^ 20 K based on the Curie-Weiss tempera¬ 
ture, Ocw = -140 K, and assuming S = 3/2). The ef¬ 
fective moment determined from Curie-Weiss analysis at 
high temperatures is 6.1 jiB] this large effective moment 
obtained at high temperatures indicates significant ther¬ 
mal population of a Jeff = 3/2 quartet at room tempera¬ 
ture. Nonetheless, the change in entropy at low tempera¬ 
tures approaches Rln2,^^ suggesting a spin-orbit coupled 
ground state Kramers doublet with Jeff = 1/2, as is often 
relevant to Co^+ materials.The degree of anisotropy 
of the effective moments in NaCaCo 2 F 7 is not yet deter¬ 
mined, but magnetization at T = 2 K and 40 K is linear 
and isotropic up to /iqH = 9 T (Appendix D).^^ 


II. EXPERIMENTAL METHOD 

We studied a 0.87 g single crystal of NaCaCo 2 F 7 using 
the MACS spectrometer at the NIST Center for Neu¬ 
tron Research.The dynamic structure factor, 5'(Q,X), 
was measured in the [HHL] reciprocal lattice plane. Two 
configurations were used for spectroscopic measurements; 
for low energy transfer scans, neutrons with a final energy 
Ef = 3.7 meV were selected, and post-sample BeO filters 
were used to remove higher harmonic contamination and 
reject neutrons for which Ef > 3.7 meV. For incident 
energies below (above) Ei = 5.2 meV, a Be filter (open 
channel) preceded the sample. For energy transfers above 
E = 8.2 meV, fixed Ef = 5.0 meV was used with post¬ 
sample Be filters and an open channel pre-sample. In 
all figures, measurements made in the various configura¬ 
tions are normalized to the same intensity units (counts 
per monitor units) using overlapping energy scans. The 
energy resolution was SE = 0.17 meV at the elastic line 
for Ef = 3.7 meV and 6E = 0.34 meV for Ef = 5.0 meV. 


III. RESULTS 

Below Tj, the elastic magnetic scattering indicates 
static short range AFM spin correlations. The subtrac¬ 
tion of 14 K from 1.7 K data reveals strong magnetic 
diffuse scattering (Fig. 1 a). There are two components 
to the this diffuse scattering; diffuse Bragg spots (not res¬ 
olution limited in Q, as will be discussed below), in addi¬ 
tion to diffuse scattering taking the shape of a “zig-zag” 
pattern underlying the peaks, i.e. the extended diffuse 
intensity is strongest along the lines connecting certain 
zone centers. This zig-zag diffuse pattern persists at fi¬ 
nite energy transfers (Fig. 4), as will be discussed below. 

The diffuse Bragg spots arise at the (111), (220), and 
(113) zone centers, but importantly not at (002) or (222) 
(see also Fig. 8 in Appendix C). These absences strongly 
constrain the frozen spin configuration. Although a short 
range structure based on collinear antiferromagnetic mo¬ 
ments might be expected based on the weak disorder 
HAFM model, the absence of the (002) diffuse Bragg spot 
rules out this scenario. The observed magnetic peaks are 
instead consistent with XY spin configurations, specifi¬ 
cally those transforming as the Fs irreducible represen¬ 
tation (IR) of the tetrahedral point group T^. The Fs IR 
admits a continuous manifold of states parameterized by 
a single angular parameter, a, which rotates the spin on 
each sublattice around its local < 111 > axis (Appendix 

A), 5,25 

At 1.7K, the elastic peak near (111) (Fig. 2 a)) can be 
described by the sum of a sharp Gaussian and a broad 
Lorentzian component. The Gaussian persists at all mea¬ 
sured temperatures; this nuclear Bragg peak is a measure 
of the instrumental Q-resolution. The Lorentzian mag¬ 
netic component gradually develops upon cooling from 
14 K (Fig. 2b), while the Full Width at Half Maximum 
(FWHM) decreases. This behavior is reminiscent of crit¬ 
ical scattering preceding a transition to an ordered state. 
From this perspective, the transition in NaGaGo 2 F 7 may 
be thought of as being preempted by freezing. Below Tj, 
the FWHM of the Lorentzian saturates at 0.12(1) A“^, 
implying a correlation length of 16.1(1) A for the short 
range magnetic order. 

Inelastic scattering is also readily observed in 
NaGaGo 2 F 7 (Fig. 3). Using the total moment sum rule 
for magnetic neutron scattering, we find that the ratio of 
elastic magnetic to total magnetic scattering is r = 0.3(1) 
in the measured region of the [HHL] plane. This is con¬ 
sistent with the ratio expected for a fully frozen (static) 
configuration arising from moments with Jeff = | (for 

which r = j ff(j^ff+i) ~ i) what is ex¬ 

pected for bare S' = | (r = 0.6). 

To interpret the inelastic magnetic neutron scattering 
we examine Q — E slices (Fig. 3 a), the energy depen¬ 
dence of constant-Q cuts (Fig. 3 b), as well as the Q- 
dependence of constant-X slices (Fig. 4). The latter 
reflect the spatial Fourier transform of spin correlations 
with a characteristic fluctuation frequency of uj = E/h. 
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FIG. 2. a) Elastic scattering (1.7 K) transverse to the (111) 
position (scan direction shown in inset). Solid line is a fit to 
a Gaussian (nuclear Bragg peak) plus Lorentzian (magnetic 
diffuse scattering). The data were acquired at fixed moni¬ 
tor counts; we report approximate counting time per setting, 
b) The temperature dependence of FWHM and area of the 
Lorentzian component of the elastic transverse scan. The fit¬ 
ting parameters for the Gaussian were fixed to their values at 
14 K. Error bars represent one standard deviation. 


From the Q — E slices of Fig. 3 a), the spectrum of 
magnetic scattering at the diffuse Bragg spots is seen to 
be gapless, while the spectrum at systematically absent 
magnetic Bragg spots (e.g. ( 002 )) is gapped and strongly 
damped. 

The constant-F^ slices reveal that the zig-zag structure 
of diffuse scattering, which is present in addition to the 
diffuse Bragg spots on the elastic line, persists to finite 
energy transfers (Fig. 4b)). For a more detailed analy¬ 
sis of these distinct spectra. Fig. 3 b) shows constant-Q 
cuts at T= 1.7 K at two locations in the [HHL] plane. 
For Q = ( 002 ) where there is no elastic magnetic peak 
(orange symbols in Fig. 3 b), the data can be fit to a 
Damped Harmonic Oscillator (DHO) spectral function 
with an excitation energy of 5.5(1) meV and a damping 
coefficient 8.7(4) meV (i.e., an overdamped mode, see 
Appendix B); we later identify this mode with collinear 
excitations out of the easy plane manifold. In contrast, 
the inelastic spectra at the diffuse Bragg positions (blue 
symbols in Fig 3b) can be fit to the sum of a Lorentzian 
relaxation function at = 0 meV with a HWHM of 
0.33(1) meV and a DHO with central energy and damp¬ 
ing fixed to the values extracted at ( 002 ) (Appendix B). 
This decomposition of the line shape of the spectrum at 
(Ill) is detailed in the right inset of Fig. 3b). 

The zig-zag pattern formed in the [HHL] plane by the 
quasi-elastic scattering (Fig.4 b)) can be associated with 
low energy states related by easy plane spin rotations 
that span the F 5 manifold. To establish this, we first 
compared the constant energy slice at 1.25 meV to the 
calculated neutron scattering intensity from a spatial av¬ 
erage of independent XY tetrahedra (Appendix C). In 
this approximation, each tetrahedron supports one choice 
from the continuous F 5 manifold. This independent XY 
tetrahedra model captures the lack of intensity near the 
( 002 ) position and the general zig-zag shape of the dif¬ 
fuse scattering. A better agreement is obtained, however, 
when collinear AFM spin components are added to each 



E (meV) 

FIG. 3. a) Magnetic scattering intensity versus energy and 
momentum transfer at T=1.7 K. The inelastic intensity near 
(002) decreases as E approaches zero (left panel), while it in¬ 
creases at the diffuse Bragg spots such as (Ill) (right panel). 
Averaging ranges: [HHO]: —0.1 < H < 0.1 (left) and [OOF]: 
— 1.1 < F < —0.9 (right), b) Guts through (Ill) and (002) 
(averaging ±0.1 r.l.u. in H and F). The solid lines are fits, 
described in the main text. The left inset shows the loca¬ 
tion of the energy cuts on a constant E = 0.25 meV map. 
The right inset shows the decomposition of the (Ill) fit into 
two components. At 1.25 meV, the same energy as the slice 
shown in Fig. 4b), both the DHO and Lorentzian components 
contribute substantially to the overall intensity (red circles). 
The data were acquired at fixed monitor counts; we report 
approximate counting time per setting. Errorbars represent 
one standard deviation. 


XY tetrahedron on the level of ^ 50% (Fig. 4 e)). This 
indicates that while the frozen state is XF-like, excita¬ 
tions for energy transfers beyond the freezing tempera¬ 
ture involve both in- and out-of plane spin components. 
In particular, the DHO mode at 5.5 meV arises from out 
of plane excitations while the quasi-elastic component is 
associated with easy plane excitations. The inset to Fig. 
3 b) shows contributions from both types of excitations 
at 1.25 meV. Beyond 4 meV, the constant energy slices 
(Fig. 4 c)) are reproduced by fully collinear AFM con¬ 
figurations on independent tetrahedra (Fig. 4 f)). An 
important shortcoming of this simple model is the in¬ 
dependent tetrahedra approximation. Specifically, the 
widths of the measured diffuse scattering for the slow 
spin fluctuations at 1.25 meV is significantly sharper in 
Q-space than predicted (Appendix C), and this is then 
evidence for inter-tetrahedron correlations. The size of 
the correlated region inferred from fitting raw data is 7.9 
A, which may be compared to the 3.65 A side length of 
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FIG. 4. Comparison of measured scattering, in counts per 
~10 seconds, at T=1.7 K (a,b,c) to simulated patterns involv¬ 
ing XY and collinear AFM static spin configurations (d,e,f). 
Measurements show diffuse scattering at a) E = 0.00 zb 0.17 
meV [14 K data subtracted, see also Fig. 1], b) = 1.25zb0.17 
meV, c) E = 4.75 zb 0.17 meV. Calculations of magnetic scat¬ 
tering from d) short range correlated cluster of XV AFM 
order with symmetry F 5 and correlation length 16 A, e) in¬ 
dependent tetrahedron model of 50% F 5 states plus 50 % 
collinear AFM spin configurations, and f) independent tetra¬ 
hedron collinear AFM model. 


been called -02 (non-coplanar) and 'ips (coplanar) (Fig. 1 
b). LRO states based on the F 5 manifold are known to 
be selected in the HAFM model upon inclusion of “indi¬ 
rect” Dzyaloshinskii-Moriya (DM) interactions,^^ despite 
an accidental ground state degeneracy admitting all val¬ 
ues of a at the mean field level. The same continuously 
degenerate manifold is also present at the mean field level 
for the XY AFM pyrochlore model,and the XT-like 
anisotropic exchange model proposed for Er 2 Ti 207 .^“^’^^ 
In all cases, the F 5 degeneracy is lifted by disorder^ and 
a LRO state is selected. The “disorder” can arise from 
thermal or quantum fluctuations, or quenched exchange 
disorder. However, in NaCaCo 2 F 7 , despite a clear mech¬ 
anism for weak bond disorder, the F 5 degeneracy is re¬ 
tained and explored by the system on short length scales 
and long time scales. In this case, exchange disorder 
does not lead to spin order, but instead to a a frozen 
spin configuration that appears to span the continuous 
F 5 manifold. 

The microscopic reason for the stabilization of XY spin 
configurations in NaCaCo 2 F 7 is not yet certain. How¬ 
ever, the spin orbit coupled Jeff = 1/2 state expected 
for Co^+ in a distorted octahedral environment could 
lead to XY anisotropy, either in the ^-tensor or the ex¬ 
change interactions, or both. The central energy of the 
damped inelastic mode (5.5 meV) may be a measure of 
the strength of the anisotropy. This should be investi¬ 
gated in the future through measurements of single-ion 
energy levels of NaCaCo 2 F 7 . Additional open questions, 
aside from quantifying the single ion anisotropy, include 
whether orbital and lattice degrees of freedom are rele¬ 
vant to NaCaCo 2 F 7 as in the related spinel compound 
GeCo204.^^ 


V. CONCLUSIONS 


the tetrahedron. Even at T = 14 K the inverse corre¬ 
lation length of the inelastic scattering is approximately 
half of that expected for independent tetrahedra. 


IV. DISCUSSION 

The above analysis of the inelastic spectrum at T = 1.7 
K suggests there are two types of dynamics in this short- 
range correlated system. The first is a relaxing process 
that is qualitatively consistent with local rotations of XY 
spin clusters through the continuous F 5 manifold with a 
relaxation time of txy = 2.02 ps. The second is a short¬ 
lived {th = 0.15 ps) inelastic mode at 5.5 meV, which is 
qualitatively consistent with collinear antiferromagnetic 
tetrahedral fluctuations, as might be expected from the 
weak disorder HAFM. 

The XY spin configurations relevant for NaCaCo 2 F 7 
at low energies are already well-studied in the pyrochlore 
literature. As shown in Appendix A, the F 5 IR can be 
decomposed into two basis vectors which have commonly 


In summary, NaCaCo 2 F 7 is the first example of a new 
class of pyrochlore single crystals based on a structurally 
ordered magnetic 3d transition metal site in a varying 
local environment created by a disordered non-magnetic 
site.^^’^^ The disordered environment leads to weak disor¬ 
der in the strong AFM interactions in NaCaCo 2 F 7 {Ocw 
= -140 K), and ultimately a low temperature freezing 
transition at Tf = 2.4 K. We have observed XY spin 
configurations forming a short range ordered state below 
Tf with a correlation length of ^ = 16 A. The low en¬ 
ergy fluctuations away from this frozen state are gapless 
to within the energy resolution of our measurement (0.17 
meV) and take on a distinctive diffuse pattern that sug¬ 
gests relaxation through a continuous manifold of local 
XY states. At higher energies, a strongly damped mode 
at 5.5 meV dominates the spectrum. The associated Q- 
dependence of the scattering intensity isconsistent with 
collinear antiferromagnetic tetrahedral fluctuations, indi¬ 
cating an XY anisotropy barrier of ^ 64 K. 

The continuous manifold of XY spin configurations 
present in NaCaCo 2 F 7 is known to collapse to an or- 
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dered state via order-by-disorder in models relevant to 
Er 2 Ti 207 as well as by DM interactions in the pyrochlore 
HAFM. However, unlike the aforementioned theoretical 
predictions, quenched exchange disorder in NaCaCo 2 F 7 
does not lead to the selection of an ordered state, but in¬ 
stead a quasi-static disordered state. Apart from the low 
energy fluctuations that appear to span the manifold, 
a prominent out of plane damped mode is observed with 
the same local structure as predicted for the Heisenberg 
model with weak exchange disorder. An intriguing aspect 
of NaCaCo 2 F 7 is the potential for ice-like correlations on 
the Na+, Ca^+ disordered sublattice. Such correlated 
disorder might be necessary to explain why NaCaCo 2 F 7 
fails to develop long range order. 
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Appendix A: Definition of states in the Ps manifold 

The sublattices of the pyrochlore lattice are described 
by the following fractional coordinates: 



The moments (pseudovectors) forming the 7/^2 and 
bases of the E 5 representation are assigned to these sub¬ 
lattices as: 



A general tetrahedral state with the symmetry of Es 
can be written as a linear combination of these sets, 

x(a) = cos a • -02 + sin a • -03 (A4) 

Assigning each “up” tetrahedron in the pyrochlore lat¬ 
tice a state x with a random value of a constitutes the in¬ 
dependent tetrahedron XY AFM state that is modeled in 
Figure 6 e). x((a) spans a continuously deformable man¬ 
ifold of XY states. These are the relevant ground states 
at the mean field level for the case of Er 2 Ti 2 07 ^“^^ or 
the HAFM model with “indirect” DM interactions.^^ 


Appendix B: Details of constant-Q lineshapes 


Figure 3 in the main text presents fits to constant- 
Q cuts, i.e. S{E). The fits to S{E) include a relaxing 
diffusive component, Sxy{E)^ and a damped harmonic 
oscillator, Sh{E). These have the well-known forms, 


Sxy{E) 


AxYE{l^n{E)) Txy 

TT + E|^ ’ 


(Bl) 


and, 

2E E 

S„(E) = A„(I + „(£)) X 

where Txy is the HWHM of the diffusing component, 
and 2Th is the damping parameter of the DHO com¬ 
ponent. The relaxation rates are then given by txy 
= 1/Txy and th = I/Th (with E’s expressed in units 
of frequency). Axy and Ah are scale factors in arbi¬ 
trary units. n{E) is the Bose-Einstein population factor, 
n{E) = {ex.p{E/kBT) — 1)“^. Ec/h is the frequency of 
the DHO mode. 



so = (l,l,2)/^/6 

§1 = (1,I,2)/V6 

5 2 = (I,1,2)/\/6 

53 = (I,I,2)/V6, 


' so = (l,I,0)/\/2 


51 = ( 1 , 1 , 0 )/V 2 

5 2 = (I,I,0)/\/2 ’ 

^ S3 = (I,1,0)/V2 

(A3) 


Axy Txy (meV) Ah Th (meV) Ec (meV) 
(002) - - 1157(52) 4.3(2) 5.5(1) 

(III) 109(2) 0.33(1) 625(17) 4.3 5.5 

TABLE 1. Parameters for fits of constant energy scans at Q 
= (002) and Q = (111) to Eqns. Bl and B2. 
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FIG. 6 . Comparison of measured diffuse scattering (a,b,c) to calculations involving XY (d,e,f), collinear AFM (j,k,l), and a 
mixture of the two types (g,h,i) of spin configurations. Diffuse scattering at a) = 0.00 zb 0.17 meV, b) E = 1.25 zb 0.17 
meV, c) E — 8.25 zb 0.34 meV. d),e),f) calculations of magnetic scattering for short range correlated cluster of F 5 XY AFM 
with correlation length 16 A (d) and independent XY tetrahedra (e,f). g),h),i) calculations involving tetrahedra with 50% XY 
and 50% collinear AFM components [g) SRO with IGA, h)i) independent tetrahedra]. j),k),l) calculations involving collinear 
AFM conhgurations. Intensities and Q-ranges in f), i) and 1) are rescaled from e), h), and k) to compare to panel c). For 
the independent tetrahedra models, each diffuse scattering pattern is computed using 50 instances of a 7x7x7 unit cell lattice 
(68600 tetrahedra in total). 


Appendix C: Supporting Neutron Scattering Data 

Here we present additional information supporting the 
conclusions from the main text. Figure 6 shows a more 
detailed comparison to three choices of models; short 
range ordered states with 16 A correlation lengths, and 
single tetrahedron states with XY or locally collinear 
character, or a mixture of these. Fig. 7 a) shows the 
widths of diffuse features as compared to the indepen¬ 
dent tetrahedra model, as well as both lower temperature 
(100 mK) and higher temperature (14 K) inelastic scans. 
The low energy diffuse inelastic scattering corresponds 
to a correlated region (^ 8 A) much larger than a single 
tetrahedron (3.5 A), at all temperatures measured, from 
100 mK to 14 K. Figure 8 a) shows an elastic scan taken 
with Ei = Ef = 1‘i.b meV at T = 1.7 K (after subtract¬ 
ing 14 K data), which reveals diffuse magnetic scattering 
throughout a larger range of Q. Note in particular the 
absence of diffuse scattering at (222). Figure 7 b) also 
shows that the inelastic scattering takes on the same pat¬ 


tern in the thermal spin liquid phase (14 K) as it does 
in the frozen phase (100 mK) (“empty can” background 
subtractions made in both panels). In Fig. 8 b and c) 
we show raw elastic scattering data (no subtraction) at 
T=1.7KandT = 14 K. 


Appendix D: Magnetization Data 

In order to investigate the possibility of an anisotropic 
^-tensor in NaCaCo 2 F 7 , magnetization measurements 
were performed at temperatures above the freezing tran¬ 
sition (T > 2.4 K), using the extraction magnetometry 
technique in a commercial physical properties measure¬ 
ment system. Measurements with the field applied along 
three different crystallographic axes were compared. The 
data taken at T = 40 K are shown in Fig. 9 a) and b). 
Only slight deviations from isotropic behavior are ob¬ 
served at the highest field strengths (^ 9 T), and these 
could easily be due to demagnetization effects for crys- 
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a) 80 



[H,H,0] [H,H,0] 

FIG. 7. a) comparison of widths of Q cuts near (220) along 
the [22L] direction, averaging over 0.28 r.l.u. in the perpen¬ 
dicular direction. The black points are taken from the dataset 
presented in Fig. 4 b) of the main text. The solid blue line is 
the width of the independent tetrahedron model represented 
in Fig. 4 e) of the main text. Green and red points are cuts 
through the 0.5 meV data shown in panel b) of this hgure. b) 
inelastic intensity oX E — 0.5 meV at 100 mK (left) and 14K 
(right). The zig-zag pattern in the 14 K (above T/) data is 
broader and more intense than at 100 mK (below T/). White 
arrow shows location and direction of cuts in a). 


tals having slightly different shapes for the different field 
orientations. Furthermore, the deviation from isotropic 
magnetization does not correspond to the expected hier- 
archy for either XY-Yike (|M(iio)| > > |M(ioo)|) 

or Ising-like (|M(ioo)| > > |M(iio)|) ^-tensors. 

For example, the magnetization of an ideal pyrochlore 
paramagnet with an XK-like ^-tensor is shown in Fig. 9 
c), using the equation. 


= gdJl^BBj{gdlJ^BJH/kBT), (Dl) 


where H is the applied magnetic field, gd is the av¬ 
erage projection of the ^-tensor onto the field direction 
(averaged over the four sites on the tetrahedron), J is 
the effective angular momentum, here taken to be 1/2 
since we may assume a spin-orbit coupled Kramers dou¬ 
blet ground state for Co^+, and Bj is the Brillouin func¬ 
tion. In Fig. 9 c) an XK-like ^-tensor was assumed, 
with gxy = 3.6 and g^ = 3.0. Although this equation is 
not expected to be valid at T = 40 K for NaCaCo 2 F 7 , 
since T < Ocw^ oiie may expect the same hierarchy of 
magnetization strengths to be observed, even in such a 
correlated paramagnetic regime. Thus, at least to within 
the demagnetization effects in these measurements, the 
^-tensor anisotropy in NaCaCo 2 F 7 is shown to be small 
on average. 
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FIG. 8. a) Elastic scattering at 1.7 K (14 K subtracted) taken with Ei = Ef = 13.5 meV. Pyrolytic graphite filters were used 
before and after the sample. Rings are due to an imperfect subtraction of the aluminum powder lines arising from the sample 
mount, b) Raw elastic data at 1.7 K (no subtraction) using Ef = 3.7 meV c) raw elastic data at 14K using Ef = 3.7 meV. 



b) 40 K 



C) T = 40K:g^y = 3.6,g^ = 3 



FIG. 9. Magnetization vs. magnetic held applied along three crystallographic directions, a), b) Magnetization data for a single 
crystal of NaGaGo 2 F 7 at T = 40 K. b) Slight differences in magnetization for the three directions can be seen a high helds, but 
the hierarchy of magnetization does not correspond to the expectations for XY-\ike or Ising-like g-tensors. c) Galculation for 
an ideal pyrochlore paramagnet (Eqn. Dl) with gxy = 3.6, Qz — 3.0, J — 1/2, and T = 40 K. 
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